rd.data.zip <- function(indata){
  
  outdata <- indata %>%
    mutate(GEOID = vlil.geoid,
           CBSA =  vlil.cbsa,
           vlil.n = round(as.numeric(as.character(vlil.fp)),3),
           lil.n = round(as.numeric(as.character(lil.fp)),3),
           mil.n = round(as.numeric(as.character(mil.fp)),3),
           hmil.n = round(as.numeric(as.character(hmil.fp)),3),
           hil.n = round(as.numeric(as.character(hil.fp)),3),
           tot.n = round(as.numeric(as.character(vlil.tot)),3),
           vli = vlil.n,
           li = lil.n - vli,
           mi = mil.n - lil.n,
           hmi = hmil.n - mil.n,
           hi = hil.n - hmil.n,
           vhi = tot.n - hil.n) %>%
    select(GEOID,
           CBSA,
           vli,
           li,
           mi,
           hmi,
           hi,
           vhi) %>%
    ## MTM added 3/1/24 
    mutate(linc = rowSums(cbind(vli,li),na.rm=T),
           minc = rowSums(cbind(mi,hmi),na.rm=T),
           hinc = rowSums(cbind(hi,vhi),na.rm=T)) %>%
    select(GEOID,
           CBSA,
           linc,
           minc,
           hinc)
  
  return(outdata)
  
}
